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Explicit  formulae  are  given  for  the  stress  discontinuities  radiated  by  a  suddenly  starting 
two-dimensional  crack  under  tension ,  in  an  ideal  elastic  body.  These  formulae  also  give , 
with  a  change  in  sign,  the  stress  discontinuities  radiated  by  a  suddenly  stopping  crack, 
and,  with  a  geometrical  deformation,  those  radiated  by  a  three-dimensional  crack.  The 
stress  in  the  primary  radiation  due  to  short  crack-jumps  is  thus  shown  to  vary  in  essentially 
the  same  manner  as  does  the  crack  speed  during  the  jump.  The  diffraction  of  the  primary 
radiation  from  one  tip  of  a  centre-crack  by  the  other  tip  produces  a  secondary  radiation, 
whose  properties  depend  mainly  on  the  nature  of  the  surface  wave  associated  with  the 
primary  radiation.  Since  this  first  diffraction  can  lead  to  crack  extension,  further  diffrac¬ 
tions  would  be  difficult  to  study  analytically.  If  it  is  possible  to  isolate  the  primary  radiation 
experimentally,  however,  it  should  give  directly  the  essential  characteristics  of  the  source. 
This  could  be  used  to  separate  source-characteristics  from  specimen  and  transducer  effects 
in  acoustic  emission  studies  of  fatigue  cracking,  or  stress-corrosion  cracking. 
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16.  ABSTRACT 

Explicit  formulae  are  given  for  the  stress  discontinuities  radiated  by  a  suddenly  starting 
two-dimensional  crack  under  tension,  in  an  ideal  elastic  body.  These  formulae  also  give, 
with  a  change  in  sign,  the  stress  discontinuities  radiated  by  a  suddenly  stopping  crack, 
and,  with  a  geometrical  deformation,  those  radiated  by  a  three-dimensional  crack.  The 
stress  in  the  primary  radiation  due  to  short  crack-jumps  is  thus  shown  to  vary  in  essentially 
the  same  manner  as  does  the  crack  speed  during  the  jump.  The  diffraction  of  the  primary 
radiation  from  one  tip  of  a  centre-crack  by  the  other  tip  produces  a  secondary  radiation, 
whose  properties  depend  mainly  on  the  nature  of  the  surface  wave  associated  with  the 
primary  radiation.  Since  this  first  diffraction  can  lead  to  crack  extension,  further  diffrac¬ 
tions  would  be  difficult  to  study  analytically.  If  it  is  possible  to  isolate  the  primary  radiation 
experimentally,  however,  it  should  give  directly  the  essential  characteristics  of  the  source. 
This  could  be  used  to  separate  source-characteristics  from  specimen  and  transducer 
effects  in  acoustic  emission  studies  of  fatigue  cracking,  or  stress-corrosion  cracking. 
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1.  INTRODUCTION 


The  technique  of  non-destructive  testing  which  relies  on  detecting  the  acoustic  emission 
(or  stress-wave  emission)  from  a  structure  under  load  does  not  yet  use  all  the  information 
contained  in  the  emission.  This  emission  is  usually  detected  by  piezo-electric  transducers  attached 
to  the  structure,  and  the  major  aim  of  current  research  is  to  derive  a  more  accurate  interpretation 
of  the  recorded  electrical  signals,  so  as  to  improve  upon  the  pragmatic  approach  which  simply 
associates  an  increasing  intensity  of  signals  with  a  deterioration  of  the  structure.  In  particular, 
this  approach  does  not  discriminate  between  the  several  sources  of  emission  which  may  be 
active,  although  their  respective  emissions  may  in  fact  show  different  characteristics,  and  thus 
be  separable  (see,  for  example,  [I]).  We  shall  be  concerned  in  this  papier  with  the  distinguishing 
features  of  the  emission  from  growing  cracks. 

The  basic  theoretical  result  is  that  a  sudden  change  in  crack  speed  gives  rise  to  stress  dis¬ 
continuities  in  the  radiated  body-waves.  Explicit  formulae  for  the  discontinuities  due  to  the 
sudden  starting  of  a  two-dimensional  crack,  under  tension  (Mode  I),  are  given  in  section  2*. 
In  section  3  we  consider  the  surface  wave  on  the  crack  faces  due  :o  the  motion  of  one  tip  of  a 
centre-crack,  and  the  secondary  radiation  due  to  the  diffraction  of  that  wave  by  the  other  crack 
tip.  A  detailed  derivation  of  the  results  given  in  sections  2  and  3  would  be  lengthy,  but  since 
it  is  not  required  for  understanding  these  results,  only  the  main  steps  are  recorded  here,  in  an 
appendix.  The  emphasis  is  placed  instead  on  the  proper  use  of  the  results  in  practice,  for  the 
theoretical  models  are,  necessarily,  very  idealized. 

In  section  4  we  show  how  Keller’s  geometrical  theory  of  diffraction  [3]  can  be  used  to 
derive  the  stress  discontinuities  radiated  by  three-dimensional  cracks,  starting  from  the  two- 
dimensional  results.  Madariaga  has  recently  performed  an  elaborate  numerical  analysis  of  the 
radiation  from  circular  shear-cracks  [4],  and  has  subsequently  shown  [5]  that  the  dominant 
features  could  be  obtained  quite  simply  from  the  relevant  two-dimensional  radiation  patterns 
(modes  II  and  III),  by  using  Keller’s  theory;  indeed,  that  approach  led  him  to  some  sharper 
insights  into  his  numerical  results.  This  makes  it  unnecessary  to  undertake  the  laborious  numerical 
calculations  for  the  analogous  tensile  cracks. 

The  theoretical  results  given  in  section  2  pertain  to  the  primary  body-wave  radiation.  It 
would  be  difficult  in  practical  applications  to  monitor  only  this  primary  radiation,  but  this  can 
be  done  in  laboratory  experiments,  if  certain  precautions  are  observed.  It  would  be  preferable 
to  capture  this  radiation  at  normal  incidence,  since  a  simple  model  can  then  be  used  for  deter¬ 
mining  the  response  of  the  piezo-electric  transducer  (section  5).  We  discuss  in  section  6  what  the 
results  of  such  an  experiment  are  likely  to  be.  Although  we  concentrate  there  on  the  longitudinal 
component  of  the  radiation,  it  should  be  noted  that  the  results  given  in  earlier  sections  are  based 
on  exact  solutions  of  the  equations  of  elasticity,  satisfying  the  proper  boundary  conditions  on 
the  crack  faces  (by  contrast  to  [6],  where  these  boundary  conditions  are  not  satisfied). 


2.  THE  RADIATION  FROM  TWO-DIMENSIONAL  CRACKS  (Mode  I) 

The  Starting  Phase 

Consider  first  the  familiar  two-dimensional  idealization  for  a  centre-cracked  plate  under 
tension  (plane  strain);  the  plate  is  assumed  to  be  of  infinite  extent;  the  crack  is  taken  to  lie  on 
the  .v-axis,  say  with  crack  tip  A  at  x  =  0,  and  tip  B  at  x  =  — /;  the  stress  at  infinity  reduces 
to  a  uni-axial  tension  a  parallel  to  the  y-axis.  The  crack  is  open,  but  in  equilibrium.  The  static 
stress  intensity  factor  is  given  by  K0  =  a(nll2y— see,  for  example,  [7].  Now  suppose  that  at 


*  Freund  [2]  discussed  briefly  the  jump  in  average  stress  across  the  longitudinal  wavefront, 
but  there  appears  to  be  an  error  in  the  numerical  results. 
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time  /  =  0  tip  A  begins  to  move  with  speed  v.  The  radiation  due  to  this  motion  will  be  called 
the  starting  phase ,  following  [5].  The  associated  wavefronts  are  shown  in  Figure  la,  where  we 
use  the  notation  P  for  the  longitudinal  wave,  5  for  the  transverse  wave,  R  for  the  surface  (Ray¬ 
leigh)  wave,  and  SP  for  the  head-wave.  The  stress  discontinuities  across  the  P-  and  5-wavefronts 
are  most  conveniently  specified  in  terms  of  polar  coordinates  r,  0.  Using  [arr]p  to  denote  the 
jump  in  the  normal  stress  a„  at  the  /‘-wavefront,  we  find  (Appendix,  §2)  that 


[<Trr]p  = 


(1  —  v)/fo 
2  (rrr)* 


V 

cl—  v  cos  0 


m, 


{( cl/ct )2  —  2  cos  0}  ( 1  +  cos  8)1 
{(cl/cr)  +  cos  6}  5  (cos  6/cl ) 


(1) 

(2) 


where  cl,  cr,  cr  are  the  longitudinal,  transverse  and  Rayleigh  wavespeeds;  v  is  Poisson’s  ratio; 
and  the  function  5  involves  an  integral  which  has  to  be  evaluated  numerically  (see  Appendix, 
equation  (A.4)).  The  following  list  of  S  (cos  8/cl),  for  values  of  cos  8  between  1  and  —1,  at 
intervals  of  0-2,  would  enable  the  reader  to  reconstruct  the  plots  in  Figure  \b,  with  the  help 
of  only  a  pocket  calculator;  0-85,  0-84,  0-83,  O' 81,  0-79,  0-76,  0-73,  0-68,  0-62,  0-51,  0-22. 
These  values,  and  all  the  numerical  results  given  subsequently,  hold  for  Poisson’s  ratio  v  =  1/3, 
which  is  a  typical  value  for  metals. 

The  important  features  of  the  result  are: 

(i)  The  strength  of  the  discontinuity  shows  a  decay  proportional  to  r~l,  which  is  charac¬ 
teristic  of  cylindrical  waves. 

(ii)  The  angular  dependence  is  best  considered  as  the  product  of  two  factors:  the  first 
leads  to  a  focussing  of  intensity  directly  ahead  of  the  crack  which  becomes  more 
pronounced  for  larger  values  of  r:  the  second  factor,  £2(8),  is  independent  of  v — it 
depends  only  on  Poisson's  ratio,  and  there  is  little  variation  for  v  in  the  usual  range 
of  interest,  0-2  v  ■  0-4.  Both  factors  are  positive  for  .8,  <  n,  with  £2  ->  0  as  8-+±n. 
Thus,  [arr]p  is  negative.  We  also  have  [a8o]p  =  {*7(1  —  *0}  [ orr\P .  The  shear  stress 
ara,  on  the  other  hand,  is  continuous  at  the  P-wavefront. 

(iii)  When  the  crack  begins  to  grow,  the  stress  intensity  factor  changes  discontinuously 
from  the  value  K0  to  the  dynamic  value  K<i  =  k(i)K0.  The  velocity  factor  k(v),  first 
derived  by  Freund  [8],  is  approximately  given  by  k(v)  =  (I  —  r /cr)  (I  —  0-95i>/cr), 
([9],  p.  802).  Thus,  [orr]p  could  have  been  related  to  the  change  in  stress  intensity 
factor  [K  ]  =  {1  —  k(v)}Kn,  rather  than  to  A0,  but  this  would  have  complicated  the 
appearance  of  equation  (I). 

At  the  5-wavefront  it  is  <rre  which  is  discontinuous,  while  the  normal  stress  components 
<7rr,  or88  are  continuous.  In  the  region  8\  <  cos ~1(—ct/cl)  —  &c,  the  jump  in  a,9  is  given  by 
(Appendix,  §2) 


farej.s  = 

(l-v)A„  «- 

4 (nrY  ct—  v  cos  0 

(3) 

/1(0)  = 

sin  20  ( ct/cl  +  cos  0)! 

(4) 

(ct/cr  +  cos  0)  5  (cos  0/cr)  ’ 

The  angular  factor  A(8)  is  positive  for  |0|  <  n/2,zero  for  |0|  =  n/2,  and  negative  for  w/2  <  |0j  <8C. 

In  the  region  [0|  >  8C,  a  head-wave  (denoted  by  SP  in  Fig.  la)  precedes  the  5-wavefront. 
It  contributes  a  logarithmic  stress-discontinuity  at  the  5-wavefront  which  may  be  derived  from 
equations  (A.  10,  II).  We  shall  not  examine  in  detail  the  behaviour  in  this  region,  because  it  is 
usual  in  laboratory  studies  to  monitor  the  acoustic  emission  ahead  of  the  crack. 

Figure  I b  shows  the  relative  magnitude  of  [<rrr]p  and  [crre].s  (for  |0i<0<-)  when  the  crack 
speed  v  =  0-5  cr,  with  the  normalization  (1  —  v)Ko /{2(Trr )J }  =  I.  The  radiation  patterns  are 
symmetrical  with  respect  to  y  =  0,  so  only  one  half  of  each  is  shown.  The  maximum  value 
of  (crre]s,  which  occurs  at  0  —  ±40°,  is  less  than  one-fourth  the  maximum  value  of  [arr]p,  which 
occurs  at  0  ===  ±85°.  The  difference  in  the  magnitude  of  these  stress  discontinuities  is  even  more 
marked  in  the  region  [0|  <30°. 
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2.2  Short  Crack-jumps 

A  stopping  phase  has  the  same  pattern  of  wavefronts,  and  the  stress  discontinuity  across 
these  has  the  same  angular  variation  as  for  a  starting  phase,  but  the  opposite  sign.  This  was 
shown  to  be  the  case  for  an  anti-plane  crack  (mode  III)  by  Eshelby  [10];  it  can  be  proved  for 
the  tensile  case  (mode  I),  as  Madariaga  [5]  has  done  for  shear  cracks  (mode  II),  on  the  basis 
of  Kostrov’s  general  solution  for  two-dimensional  crack  problems  [II],  Thus,  equations  (1) 
and  (3),  without  the  first  minus  sign,  also  hold  for  a  stopping  phase,  but  v  now  denotes  the 
crack’s  speed  just  before  it  stops,  and  Ko  must  be  replaced  by  a  factor  K*  which  depends  on 
the  precise  history  of  the  crack’s  motion  (cf.  [12],  §2).  In  practice  one  is  often  concerned  with 
crack  jumps  which  are  much  shorter  than  the  current  crack  length,  and  it  is  then  a  reasonable 
approximation  to  take  K*  equal  to  the  static  stress  intensity  factor  for  the  current  crack  length. 

There  is  no  need  to  consider  in  detail  what  happens  behind  the  wavefronts  if  one  is  only 
interested  in  short  crack-jumps:  for  practical  purposes,  the  variation  of  the  compressive  stress 
in  the  P-pulse  (and  of  the  shear  stress  in  the  5-pulse,  ahead  of  the  crack)  is  approximately  the 
same  as  the  variation  in  crack  speed  during  the  jump.  For  example,  if  the  crack  both  starts 
and  stops  abruptly,  the  radiation  ahead  of  the  crack  will  consist  of  two  approximately  rectangular 
stress-pulses,  whose  duration  will  depend  on  the  length  of  the  jump,  and  whose  amplitude  will 
depend  on  /Co  and  the  crack  speed  v  during  the  jump.  It  can  be  expected  that  in  practice  the 
pulses  will  not  be  that  sharp,  for  at  least  two  reasons:  first,  changes  in  crack  speed  are  more 
likely  to  be  only  “effectively  discontinuous”,  relative  to  an  appropriate  time-scale,  rather  than 
strictly  discontinuous;  and  secondly,  in  the  two-dimensional  idealization  all  points  on  a  straight 
crack-front  move  together,  whereas  in  practice  it  is  more  likely  that  different  portions  of  the 
crack  front  move  at  different  speeds,  and  even  at  different  times. 

In  spite  of  these  smoothing  effects,  the  body-wave  radiation  from  growing  cracks  can  still 
be  expected  to  show  “effective  discontinuities”  in  stress,  due  to  “relatively  sudden”  variations 
in  crack  speed,  which  occur  typically  when  a  crack  begins  to  grow,  or  when  a  growing  crack  stops. 


2.3  Correlation  with  the  Stress  Intensity  Factor 

The  theoretical  results  for  moving  cracks  show  that  the  stress  variation  across  a  radiated 
wavefront  is  essentially  the  same  as  the  variation  of  the  stress  intensity  factor  at  the  time  of 
generation  of  the  wave  (cf.  section  2. 1  (iii)).  This  is  also  true  for  stationary  cracks.  For  example, 
when  a  tensile  load  is  suddenly  applied  to  a  centre-cracked  plate,  the  stress  intensity  factor  K 
increases  in  direct  proportion  with  t \  if  the  loading  step-pulse  strikes  the  crack  at  normal  incidence 
at  time  t  =  0.(K continues  to  increase  until  the  diffracted  wave  from  either  tip  reaches  the  other, 
when  it  begins  to  decrease,  and  then  to  oscillate  about  its  eventual  static  value,  [13].)  It  can  be 
shown  that  now  the  radiated  (or  “diffracted")  wavefronts,  from  the  crack  tip  at  .v  =  0,  carry 
a  change  in  stress  which  is  initially  proportional  to  t!,  if  r  denotes  the  time  elapsed  since  the 
arrival  of  the  P-  or  5-wavefront  at  any  given  point  (r,  0)  ahead  of  the  crack  (tf!  <w/2,  r=(r— r/ci) 
or  ( t—rjcr )  respectively).  It  is  sufficient  in  fact  to  verify  that  this  is  the  case  for  a  semi-infinite 
crack,  by  starting  from  the  results  given  in  [14],  §9.6.3. 

Notice  that  in  this  case,  with  a  stationary  crack  tip,  K  does  not  show  a  step-discontinuity, 
even  though  the  incident  pulse  is  discontinuous.  It  follows  that  if  the  starting  phase,  or  the 
stopping  phase,  from  one  crack  is  diffracted  by  another  neighbouring  crack,  the  stress  across 
the  diffracted  wavefronts  will  be  continuous.  Thus,  stress  discontinuities  (or  “effective  discon¬ 
tinuities”)  are  distinguishing  features  of  the  radiation  from  growing  cracks. 


3.  THE  DIFFRACTION  OF  THE  STARTING  PHASE 

Figure  2  shows  the  displacement  of  the  lower  half  of  the  crack  surface  when  the  P-wavefront 
of  the  starting  phase  from  crack  tip  A  has  just  reached  tip  B,  which  we  suppose  to  have  remained 
stationary.  The  only  prominent  feature  is  a  small  hump  travelling  to  the  left  with  speed  cr ; 
even  that  is  discernible  only  because  we  have  assumed  that  tip  A  moves  at  the  relatively  high 
speed  p  =  0-5  a.  (see  Appendix,  §3). 

The  question  now  is  to  determine  the  variation  of  the  stress  intensity  factor  at  B  due 
to  the  diffraction  of  the  starting  phase.  This  will  give  as  well  the  stress  variation  across  the 
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diffracted  wavefront  radiated  from  B,  according  to  the  general  result  explained  in  section  2.3. 
The  answer  applies  also  to  the  case  where  both  tips  A  and  B  start  to  move  and  stop  shortly 
after,  before  the  starting  phase  from  either  tip  has  reached  the  other.  We  find  that 

*/ 

AKn{t)  ------  -Ko  g(u)/(llcn  -  u)  du,  (5) 

J 

where  g(«)  is  a  regular  positive  function  (Appendix,  §4).  Thus,  AKr  -  0  until  t  —  Ijcu  when 
the  /’-wavefront  from  A  reaches  B\  it  then  takes  a  negative  value,  of  gradually  increasing  magni¬ 
tude,  and  tends  logarithmically  to  a  negatively-infinite  value  as  t  approaches  //c«,  the  arrival 
time  of  the  surface  wave  associated  with  the  starting  phase.  This  singular  behaviour  is  due  to 
the  existence  of  a  surface  wave.  Quite  a  different  behaviour  would  be  predicted  for  the  corres¬ 
ponding  anti-plane  problem,  where  there  is  only  one  type  of  body-wave  and  no  surface  wave. 
There,  the  theoretical  result  is  that  AKr  =  0  until  t  —  IJct.  and  then  it  increases  smoothly 
and  monotonically— see  (A.  1 9). 

With  the  tensile  crack,  the  stress  intensity  factor  at  B  falls  to  zero  when  AKr  =  —  Ko 
(i.e.  before  t  =  I/cr).  The  crack  is  then  effectively  closed  in  the  neighbourhood  of  B.  The  theore¬ 
tical  result  is  not  valid  beyond  that  time,  for  it  would  predict  that  the  crack  faces  press  against 
one  another,  whereas  in  its  derivation  it  is  supposed  that  this  does  not  happen.  However,  this 
difficulty  does  not  arise  for  the  diffraction  of  a  stopping  phase,  and  the  theoretical  prediction 
then  is  that  A K  will  increase  logarithmically  to  a  positively-infinite  value.  According  to  the  theory 
of  quasi-brittle  fracture  (see  [7]),  the  crack  will  begin  to  extend  at  tip  B  if  Kb  exceeds  the 
fracture  toughness.  It  may  stop  shortly  after,  if  the  driving  force  is  then  below  the  critical  value. 
Thus,  a  stopping  phase  from  one  tip  can  lead  to  a  sequence  of  secondary  radiation,  due  to  the 
diffraction  of  successive  stopping  phases.  This  would  be  difficult  to  keep  track  of  analytically. 

The  theoretical  results,  however,  reflect  the  idealizations  of  the  model.  In  addition  to  those 
discussed  in  section  2.2,  there  is  the  idealization  that  the  crack  faces  are  perfectly  smooth:  the 
surface  wave  then  propagates  unattenuated.  In  practice,  the  crack  faces  are  rough,  and  will  scatter 
the  high-frequency  component  of  the  surface  wave,  converting  it  into  body-wave  radiation.  We 
note  also  that  the  logarithmic  approach  to  infinity  is,  in  mathematical  terms,  less  singular  than 
the  “inverse  square-root"  singularity  reported  by  Freund  [15]  (see  also  [16],  equation  (35))  for 
the  case  where  point-forces  are  suddenly  applied  to  the  faces  of  a  crack.  Freund  [15]  showed 
that  the  singularity  is  subdued,  and  virtually  disappears,  if  the  forces  are  applied  gradually.  One 
can  expect  a  similar  smoothing  to  occur  if  the  crack  accelerates  gradually  to  the  speed  r,  but 
the  result  in  the  present  case  cannot  be  obtained  by  a  simple  convolution:  it  requires  a  detailed 
tracking  of  the  location  of  the  crack  tip  at  various  retarded  times  (cf.  [12],  §4,  for  a  similar 
anti-plane  problem). 

Although  it  is  not  possible,  on  purely  theoretical  grounds,  to  determine  precisely  what  will 
happen  in  practice,  it  should  be  emphasized  that  the  surface  wave  associated  with  the  starting 
phase  is  not  a  theoretical  artifact,  but  a  real  phenomenon.  In  particular,  it  is  noted  by  Shmuely 
et  al.  [17]  as  being  responsible  for  the  seemingly  premature  crack-arrest  observed  for  certain 
geometrical  configurations,  in  their  experiments. 

4.  THE  RADIATION  FROM  A  CIRCULAR  CRACK 

Suppose  that  at  time  /  0  a  circular  crack  appears  on  the  plane  ;  =  0,  in  a  field  of  uniform 

tension  parallel  to  the  e-axis.  The  crack  grows  with  constant  speed  r  from  zero-radius  to  a 
radius  a,  and  then  suddenly  stops.  For  the  purposes  of  illustration,  it  will  be  sufficient  to  consider 
only  the  longitudinal  component  of  the  radiation.  The  variation  of  stress  immediately  behind 
the  spherical  /’-wavefront  generated  by  the  appearance  of  the  crack  is  the  same  as  it  would  be 
if  the  crack  did  not  stop  but  continued  to  grow  with  speed  r.  It  can  be  verified,  from  the  known 
solution  for  the  self-similar  circular  crack  (see.  for  example.  [18]),  that  the  normal  stress  com¬ 
ponents  are  proportional  to  r,  the  time  elapsed  since  the  arrival  of  the  wavefront. 

The  dominant  feature  of  the  radiation,  however,  is  the  stress  discontinuity  carried  by  the 
stopping  phase,  as  it  is  clearly  shown  by  Madariaga's  work  for  the  analogous  shear  crack  [5). 
The  manner  in  which  this  discontinuity  is  registered  at  any  point  may  be  derived  from  Keller's 
geometrical  theory  of  diffraction  [3]  as  follows.  We  imagine  that  each  infinitesimal  element  of 
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the  crack  front  generates,  on  stopping,  a  fan  of  rays  in  a  plane  normal  to  that  element.  With 
a  circular  front,  only  two  of  these  rays  pass  through  any  given  field-point  Q:  they  come  from 
that  point  on  the  front  which  is  closest,  and  that  which  is  furthest,  from  Q  (Fig.  3).  The  ray 
from  the  closest  point  carries  the  step-discontinuity 


(I  -  v)  K  _ 

2|»K1(I  l-  RiIpiW 


CL  — 


V 

v  cos  (9 
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which  is  derived  from  the  two-dimensional  result  given  in  equation  (I)  by  making  the  following 
substitutions,  using  the  notation  shown  in  Figure  3: 

(i)  Ko  is  replaced  by  K,  which  is  approximately  equal  to  the  static  stress  intensity  factor 
for  a  circular  crack  of  radius  a  (cf.  [5],  §8); 

(ii)  6  is  replaced  by  /S; 

(iii)  r  s  is  replaced  by  j/?i(l  f  Rilpi)}  !,  with  pi  —  ajcos  fi. 

The  last  change  is  the  most  important  one,  and  it  amounts  to  a  geometrical  deformation 
to  account  for  the  curvature  of  the  crack  front.  For  large  R.  the  amplitude  of  the  discontinuity 
decays  as  R  l,  as  it  would  for  a  spherical  wave.  The  difference  in  sign  between  equations  (I) 
and  (6)  is  due  to  the  fact  that  we  are  now  considering  a  stopping  phase.  The  stress  component 
denoted  by  <r  in  (6)  is  the  normal  stress  parallel  to  the  ray;  the  discontinuity  in  the  normal  stress 
at  right-angles  to  the  ray  is  also  given  by  (6)  but  with  ( 1  —  v)  replaced  by  v  (cf.  the  relation 
between  aT,  and  aee  in  section  2.1).  Since  Keller's  theory  is  valid  asymptotically  for  high  frequ¬ 
encies,  it  can  be  expected  to  give  accurately  the  values  of  these  step-discontinuities,  and  it  is 
evident  from  Madariaga’s  work  that  this  is  indeed  the  case.  The  discontinuity  carried  by  the  ray 
from  the  furthest  point  is  derived  in  a  similar  manner,  but  with  6  replaced  by  (n  —  y).  The 
delay  between  these  two  rays  gives  a  useful  measure  of  the  effective  width  of  the  /*-pulse  of  the 
primary  radiation  (except  for  a  =  w/2,  where  there  is  a  caustic).  In  the  far  field,  i.e.  for  Rp  a, 
this  delay  is  approximately  equal  to  (2a  cos  a/cc). 

The  diffraction  of  the  stopping  phase  will  give  rise  to  a  secondary  radiation,  whose  amplitude 
can  again  be  derived  on  the  basis  of  Keller’s  theory  (because  the  symmetry  of  the  present  model 
implies  that  the  energy  redistributes  itself  uniformly,  after  the  rays  cross  the  caustic  at  the  origin), 
but  this  will  lead  to  the  same  problems  of  interpretation  as  were  dealt  with  in  the  previous 
section.  Also,  the  detailed  calculations  for  the  transverse  component  of  the  primary  radiation 
will  be  complicated  in  the  region  where  the  head-wave  precedes  the  S-wavefront. 

The  suddenly  stopping  circular  crack  may  be  considered  as  an  idealized  model  for  the 
fracture  of  brittle  inclusions.  Its  most  serious  idealization  is  that  the  whole  crack  front  stops 
simultaneously  and  suddenly.  In  practice,  it  is  more  likely  that  the  fracture  initiates  near  the 
inclusion/matrix  interface,  so  that  different  portions  of  the  crack  front  come  to  rest  at  different 
times.  The  implication  here,  as  in  sections  2  and  3,  is  that  the  radiation  from  the  model  has 
sharper  features  (viz.  stress  discontinuities),  and  a  larger  amplitude,  than  what  is  likely  to  be 
observed  in  practice.  Nevertheless,  the  theoretical  results  can  provide  a  useful  guide  in  designing 
experiments  and  interpreting  their  results  (see  section  6). 


5.  THE  RESPONSE  OF  A  PIEZO-ELECTRIC  TRANSDUCER  TO  A  NORMALLY 
INCIDENT  P-PULSE 

There  is  a  simple  one-dimensional  model  for  determining  the  response  of  a  piezo-electric 
transducer  to  a  normally  incident,  plane  /’-pulse.  Although  the  model  has  been  known  for  some 
time  (see  [19]),  its  potential  usefulness  appears  to  have  been  overlooked.  It  consists  of  an  elastic 
slab,  of  uniform  thickness  d,  sandwiched  between  two  elastic  half-spaces.  The  slab  gives  rise 
to  a  voltage  proportional  to  the  relative  normal  displacement  of  its  faces;  thus,  the  voltage  is 
proportional  to  the  integral  of  the  normal  strain ,  and  also  to  the  integral  of  the  normal  stress , 
across  the  thickness  of  the  slab.  The  essential  parameters  are  the  acoustic  impedances,  which 
will  be  denoted  by  Ho  for  the  slab,  W\  for  the  half-space  in  which  the  pulse  originates,  and  Wz 
for  the  backing  (to  conform  with  the  notation  of  [19],  p.  1 16).  In  applying  the  model  to  a  parti¬ 
cular  transducer,  Ho  is  taken  to  be  equal  to  pncn,  where  po  is  the  density  of  the  piezo-electric 
material,  and  <n  is  the  speed  of  longitudinal  waves  across  the  thickness;  the  value  for  PZT  5A, 
which  may  be  considered  as  typical,  is  Ho  —  34  (kg  m  3  x  km  s  ')*. 

*  I  am  indebted  to  Mr  I.  G.  Scott  for  supplying  this  information. 
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W\  is  the  acoustic  impedance  of  the  specimen;  for  example,  W\  =  17  (kg  nr*  x  kms1)  for 
aluminium.  It  is  usually  the  case  that  the  impedance  Wi  of  the  backing  material  is  less  than  Wi, 
so  that  there  is  reflection  with  a  change  of  sign  at  the  back  interface.  This  leads  to  the  charac¬ 
teristic  ring-down  of  the  recorded  voltage  in  response  to  short  pulses  (see  [19],  Fig.  1.118,  in 
which  Wo/  W2  =  9). 

To  use  the  model  we  need  its  impulse  response  /(/).  This  is  not  given  in  [19],  but  it  is  not 
difficult  to  show  that  it  can  be  represented  analytically  by  the  following  convolution  product 
(denoted  by  an  asterisk): 


/(/)  = 


(Ri  R2)n  8  [(c0t/2d)  -  «)] 


T  =  2W0KW0  + 


0) 


(I,  0  <  t  <  d/co , 

Ri,  d/co  <  I  <  2d  I  Co, 

0,  t  <  0,  /  >  2d/co , 

/?1  =  {Wi  -  W0)/(W0  +  Wl),R2  =  (W2-  W0)!(Wo  +  I V2). 

8  denotes  the  Dirac  delta-function,  k  denotes  the  dynamic  piezo-electric  constant;  its  actual 
value  is  not  required  for  reconstructing  the  shape  of  the  incident  pulse.  The  step  response  A  (t), 
which  is  the  integral  of  1(1),  is  an  oscillating  function  which  approaches  asymptotically  the 
limiting  value  kT (\  +  /?2>/(  1  —  R1R2),  as  shown  in  Figure  4.  The  response  to  a  rectangular 
pulse  of  unit  amplitude  and  duration  At  is  A  (t)  —  A  (t  —  At). 

The  model  can  be  expected  to  give  accurately  only  the  initial  response  of  a  transducer, 
because  the  finite  dimensions  of  the  transducer,  and  the  cross-coupling  in  piezo-electric  materials, 
will  play  an  increasingly  important  role  in  determining  the  response.  Conversely,  given  an  actual 
transducer  output,  one  can  expect  to  reconstruct  accurately  only  the  beginning  of  the  stress 
pulse.  But  this  is  often  sufficient  to  identify  the  source,  as  noted  in  [1],  p.  98.  Also,  if  the  duration 
of  the  pulse  is  not  much  greater  than  the  transit  time  d/c0  (which  is  typically  0-5  ps),  one  can 
expect  to  reconstruct  the  whole  pulse  accurately,  and  for  this  it  is  sufficient  to  retain  only  the 
first  few  terms  of  the  series  in  equation  (7).  Figure  4  only  required  the  first  four  terms;  this 
should  be  sufficient  in  practice.  The  reconstruction  involves  a  de-convolution  which  generally 
has  to  be  performed  numerically.  This  may  be  done  directly  in  the  time-domain  (as  in  [1]),  or 
in  the  frequency-domain,  where  de-convolution  is  equivalent  to  a  division,  by  using  the  recently 
developed  Fast  Fourier  Transform  algorithm  [20]. 

The  experimental  results  shown  in  Figure  15  of  [I]  suggest  that  the  model  can  in  fact  be 
useful  for  stress  pulses  of  much  longer  duration  than  d/ca.  That  figure  shows  the  surface  dis¬ 
placement  as  measured  by  a  capacitive  transducer,  and  the  response  of  a  piezo-electric  transducer, 
due  to  a  normally  incident  pulse.  The  overall  shape  of  the  displacement/time  record  is  close 
to  that  of  a  ramp-function,  with  a  rise-time  of  30  ps\  if  it  were  precisely  a  ramp-function,  the 
corresponding  stress  pulse  would  be  rectangular  and  of  30  ps  duration.  The  response  of  the 
piezo-electric  transducer  is  in  fact  similar  to  that  shown  in  our  Figure  4;  indeed,  one  can  account 
for  the  finer  details  of  the  initial  response  by  taking  into  account  more  precisely  the  actual  stress 
variation. 

It  is  important  to  note  that  the  model  assumes  perfect  acoustic  contact  between  the  trans¬ 
ducer  and  the  specimen.  This  is  usually  ensured  in  practice  by  applying  a  slight  pressure  to  the 
transducer  in  attaching  it  to  the  surface.  Also,  it  is  not  difficult  to  allow  for  the  presence  of  a 
viscous  couplant  between  the  transducer  and  the  specimen,  within  the  framework  of  a  one¬ 
dimensional  model.  The  essential  point  is  that  while  a  proper  modelling  of  the  long-term  response 
of  a  piezo-electric  transducer  would  be  complicated,  its  initial  response  to  a  normally  incident 
pulse  can  be  determined  satisfactorily  from  a  simple  model.  This  can  be  used  to  advantage  in 
designing  experiments,  as  illustrated  in  the  next  section. 


6.  DISCUSSION 

To  simplify  the  discussion  we  shall  consider  a  particular  specimen,  shown  in  Figure  5a,  and 


6 


concentrate  on  the  emission  from  one  increment  of  crack  growth.  For  example,  with  a  fatigue 
crack,  it  would  be  convenient  to  isolate  the  emission  produced  by  an  increase  in  the  maximum 
load  for  one  cycle,  after  a  steady-state  has  been  reached  under  constant  amplitude  loading. 

(There  are  differences  of  opinion  on  whether  the  acoustic  emission  from  fatigue  cracks  is 
due  mainly  to  crack  growth  or  to  dislocation  sources  (cf.  [21,  22]),  or  even  whether  the  crack 
growth  occurs  near  the  peak  of  the  loading  half  of  the  cycle  or  during  unloading  (cf.  [23,  24]). 
The  point  of  view  adopted  here  may  be  summarized  as  follows.  Consider  an  actual  compact 
tension  specimen,  and  an  ideally  elastic  model  of  it.  Let  the  crack  length  in  the  model  be  adjusted 
according  to  the  plastic-zone  correction  for  small-scale  yielding  (see  [7],  §3.12;  [25]).  Then, 
during  the  early  part  of  the  overload  cycle,  the  elastic  stress  field  beyond  a  small  region  surround¬ 
ing  the  crack  tip  will  be  the  same  in  the  specimen  as  in  the  model.  To  maintain  the  correspondence 
as  the  peak-load  is  approached,  one  must  extend  the  crack  in  the  model.  What  is  not  known 
in  advance  is  the  variation  in  crack  speed  during  this  extension  which  will  give  the  best  corres¬ 
pondence  between  the  theoretical  and  the  actual  radiation.  Once  that  is  determined  experi¬ 
mentally,  one  has  a  macroscopic  characterization  of  the  source  which  does  not  call  into  question 
the  microscopic  mechanisms  of  deformation  and  rupture  which  are  responsible  for  the  actual 
emission.) 

In  actual  experiments,  care  should  be  taken  to  absorb  as  much  as  possible  of  the  primary 
radiation  when  it  reaches  the  edge  of  the  specimen,  especially  the  surface  wave  discussed  in 
section  3  when  it  reaches  the  corners  C  in  Figure  5a,  at  least  until  the  primary  radiation  is  well 
understood.  We  suppose  that  the  experiment  is  set  up  to  record  this  primary  radiation,  say 
at  an  angle  of  30°  to  the  direction  of  the  crack.  The  specimen  in  Figure  5a  is  a  compact  tension 
specimen  (like  that  used  in  [22]),  but  with  one  corner  cut  off  so  that  the  transducer  at  location  A 
receives  the  radiation  at  normal  incidence.  The  specimen  is  supposed  to  be  sufficiently  thick  to 
accommodate  fully  that  transducer. 

The  normal  stress  arr  carried  by  the  P-pulse  can  now  be  reconstructed  from  the  response 
of  a  piezo-electric  transducer  at  A ,  by  using  the  model  of  section  5,  or  the  corresponding  normal 
surface  displacement  at  A  can  be  monitored  by  using  a  capacitive  transducer,  as  in  [I],  In  either 
case,  the  properties  of  the  /’-pulse  should  be  sufficient  to  characterize  the  source,  for  short 
crack-jumps  (cf.  section  2.2). 

The  transducer  at  B  is  placed,  as  usual,  on  one  of  the  broad  faces  of  the  specimen,  say 
the  face  lying  in  the  plane  z  =--  h .  the  other  face  lying  in  the  plane  r  =  —h.  It  would  be  difficult 
to  predict  the  response  at  B,  and  what  follows  is  only  a  qualitative  discussion  of  what  can  be 
expected  there  if  the  P-pulse  registered  at  A  is  a  narrow  rectangular  pulse,  approaching  an 
impulse.  The  first  difficulty  is  that  the  results  of  section  2  are  valid  only  if  a  state  of  plane  strain 
prevails  throughout  the  specimen,  but  this  is  not  the  case  on  the  faces  c  =  ±h.  The  most 
important  modification  which  must  be  made  is  to  cancel  the  stress  ozz  (=  v(oTT  -r  tree)),  pre¬ 
dicted  by  the  plane-strain  solution,  over  an  expanding  quasi-annular  region,  representing  the 
intersection  of  the  two-dimensional  circular  wavefronts  of  the  P-pulse  with  the  surfaces  s  =  ±h. 
The  expected  normal  displacement  uz  on  :  =  h  shown  in  Figure  5b  was  obtained  by  smoothing 
the  singularities  of  the  theoretical  solution  for  uz  due  to  an  expanding  circular  ring-load  on  a 
half-space  [26],  but  retaining  a  sense  of  the  relative  magnitudes  of  (he  various  components 
of  this  surface  “bulge".  It  is  not  called  a  surface  wave  because  the  leading  edge  clearly  travels 
with  the  speed  of  longitudinal  waves,  which,  in  metals,  is  approximately  twice  the  speed  cr  of 
surface  waves.  On  the  other  hand,  it  is  also  clear  from  Figure  5b  that  the  major  hump  does 
travel  with  speed  cr,  and  it  can  properly  be  called  a  surface  wave:  in  particular,  there  would  be 
the  characteristic  retrograde-elliptical  particle  motion  associated  with  that  hump  (cf.  [14],  §5.1 1). 
The  response  of  a  piezo-electric  transducer  at  B  will  depend  critically  on  the  manner  in  which 
it  is  attached  to  the  surface,  for  example  on  whether  it  is  perfectly  bonded  to  the  surface,  or 
kept  at  a  certain  distance  from  it  due  to  the  use  of  a  couplant.  Whatever  the  precise  conditions 
may  be,  we  can  anticipate  that: 

(i)  the  apparent  duration  of  the  pulse  will  depend  on  the  size  of  the  transducer,  being 
longer  for  a  larger  transducer; 

(ii)  unlike  the  /’-pulse,  the  surface  bulge  will  not  be  narrow,  and  indeed  its  width  will 
increase  with  distance  from  the  source,  since  the  leading  edge  travels  faster  than  the 
surface  wave  component; 

(iii)  in  using  the  theoretical  solution  for  a  ring-load  on  a  half-space,  we  have  ignored  the 
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reflections  back  and  forth  between  the  faces  z  =  ±h,  which  will  further  contribute  to 
the  apparent  broadening  of  the  pulse;  and 

(iv)  if  B  is  sufficiently  far  from  the  crack  tip,  the  body-waves  generated  by  cancelling  azz 
on  z  =  —h  will  reach  B  before  the  surface  wave  travelling  on  z  =  h. 

These  comments  also  apply  to  the  response  of  a  capacitive  transducer  at  B. 

The  crux  of  the  discussion  is  this:  with  the  transducer  at  A  we  are  dealing  with  an  essentially 
two-dimensional  situation;  with  the  transducer  at  B  we  are  inextricably  in  a  difficult  three-dimen¬ 
sional  situation.  Although  in  the  practical  application  of  acoustic  emission  monitoring  the 
transducer  will  generally  be  in  location  B ,  in  laboratory  experiments  an  additional  transducer 
can  be  used  in  location  A.  This  should  allow  one  to  reconstruct  the  /‘-pulse  of  the  primary 
body-wave  radiation,  which  should  characterize  the  source  and  thus  provide  a  secure  basis  for 
the  interpretation  of  the  signals  registered  at  B.  In  particular,  this  could  be  used  to  separate 
source-characteristics  from  specimen  and  transducer  effects  in  acoustic  emission  studies  of  fatigue 
cracking,  or  stress-corrosion  cracking. 
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APPENDIX 


1.  Introduction.  Notation 

The  problems  to  be  considered  below  are  solved  by  integral  transform  methods:  the  basic 
techniques  are  explained  in  [14],  chapters  7  and  9;  the  major  innovation  required  here  is  due  to 
Kostrov  [II],  We  shall  adopt  the  notation  used  in  our  previous  work  [16],  subsequently  referred 
to  as  I.  In  particular,  a  circumflex  is  used  to  denote  a  Laplace  transform  over  time,  with  para¬ 
meter  s,  and  the  corresponding  capital  letter  is  used  when  the  transformation  over  time  is  followed 
by  a  bilateral  Laplace  transformation  over  x,  with  parameter  s\-  a,  b,  c,  d  denote  the  reciprocals 
\/cL,  \/ct,  1  Icr,  1  /*?;  A,  n  the  Lame  constants.  The  crack  problems  may  be  reduced  to  mixed 
boundary-value  problems  for  a  half-space  (y^O),  and  the  condition  of  symmetry  (crxv(>'=0)=0) 
leads  to  the  relation  (I:  (7)) 

Uy(x,  O,  s)  —  —K  {A/(x)/j}  2yy(.X’  0,s),  ( A.  1 ) 

X  =  (1  —  v)//x  ,  (A.2) 

on  the  boundary  ( y  =  0)  of  the  half-space.  In  the  solution  of  crack  problems  it  is  necessary  to 
factorize  the  kernel  M  (\)  as  follows  (I,  §4): 

M  =  M~(x)  =  M+(~x) . 

MHx)  =  (a  +  x)‘/{(c  +  x)S  (*)} ,  (A. 3) 

S  (x)  =  exp 

A  supercript  +  attached  to  a  function  of  x  indicates  that  the  function  is  “analytic  to  the  right”, 
and  consequently  its  inverse  is  identically  zero  for  x  <  0. 

2.  The  starting  phase 

The  appropriate  boundary  conditions  are: 

Vyy{o  <  X  <  Vt,  O,  t  ^  o)  =  -Kolibrxy  , 

Uy„(x  <O,O,t)  =  0, 

Uy(X  >  vt,  o,  t)  =  0  . 

Notice  that  we  are  using  only  the  singular  part  of  the  prior  static  field  ahead  of  the  crack:  this 
simplifies  considerably  the  expression  for  the  dynamic  stress  field  generated  by  the  motion 
(cf.  the  discussion  of  the  static  factor  K *  in  [12],  §2),  but  it  gives  accurately  the  stress  step- 
discontinuities  at  the  wavefronts. 

The  main  steps  are  the  following: 

(i)  First,  the  double  transform  of  Kostrov’s  pseudo-stress  P(\,  s),  defined  in  I:  (52),  is 
derived.  This  proves  to  be  the  same  as  the  double  transform  of  the  actual  stress  for 
the  corresponding  anti-plane  problem  (cf.  [5],  §4). 

(ii)  From  this  we  find  the  double  transform  of  the  actual  stress  on  y  =  0  to  be 

Zyt Ax,  o,  s)  =  P  (x,  s)/{{b  +  x)1  M<(x)),  (by  definition), 

=  —  Kn/{\/2  v3  \d  +  x)  •  (A.5) 

(iii)  Using  the  equations 

s2  4>  (x,  o,  s)  =  -(*c/62)  (y/«)  •*) .  (A.6) 


- I  arctan 


4m2  (u2  —  a2)1  (b2  —  u2y  du 


(2m2—  b2)2 


u  +  X 


(A.4) 


V  (x>  s)  =  -(2  x«/y)  *  (x.  o>  s) .  (A. 7) 

(where  a  denotes  a  branch  of  (a2  —  x2)1  with  positive  real  part,  y  =  (2*2  —  b2),  cf.  I : 
(6)),  we  can  derive  the  double  transforms  of  the  displacement  potentials  <j>,  i/>  and 
proceed  to  perform  their  inversion  by  the  Cagniard — de  Hoop  technique  ([14],  §7.9). 
With  the  notation 


£  =  — (f/r)  cos  8  +  /  {(///•  )2—  a2}1  sin  8 , 

where  £  is  the  value  of  x  along  the  Cagniard  contour  in  the  upper  half-plane,  and 
(r,  8)  are  the  polar  coordinates  of  the  point  (x,  y),  we  have 


£,,  =  /'(a2-  £2)‘/(/2—  a2  r2)1, 


l*x.i],-c  =  U2-  a2  r2)1  Re  {«*  , 


so  that  the  inversion  integral  with  respect  to  s\  is  brought  into  the  form 


<t>  ( r ,  8,  s)  =  s~612  (A/n) 


Hr,  8,  t) 

(,2_  a2  r2)i 


e ~ ,l  dt , 


(A. 8) 


where  A  —  ( x/b 2)  (A'o/v/2),  and/(r,  8,  t)  is  the  value  of  Re  [s5'2  (x<P/A ]  along  the  Cagniard 
contour  x  =  £•  A  similar  expression  can  be  derived  for  0  (r,  8,  s)  in  the  region  |0|  <8C. 


(iv)  The  inversion  with  respect  to  s  now  proceeds  by  using  the  convolution  theorem.  It  is 
more  convenient  to  carry  out  this  inversion  for  the  time  derivatives  because 

these  show  directly  the  stress  variation  at  the  wavefront  (see  below).  We  find  that 


<£,u  (r,  8,  t)  =  Ait  3/2 


'  fir,  8 ,  u)  du 

ar  ("2  -  a2  r2Y  (r  -  «)*  ’ 


so  that,  by  virtue  of  the  result 


(«(r  —  «)}-*  du  =  v 


which  holds  for  all  r,  and  in  particular  for  r  0+,  is  found  to  have  a  step-dis¬ 
continuity  at  the  /’-wavefront,  given  by 


=  M/(2war)»}/(r,  8,  ar). 


(A.9) 


It  will  be  noticed  that /(r,  8,  ar),  which  gives  the  angular  variation  of  the  discontinuity, 
is  in  fact  equal  to 


Re  [s5/2  a.<l>l A  ]  =  Re  [yM-(x)l(d  +  x)]. 

evaluated  at  the  point  x  =  a  cos  8  where  the  Cagniard  contour  intersects  the  real 
axis.  This  angular  variation  is  the  same  as  that  given  in  equations  (I)  and  (2)  of  the 
main  text. 


(v)  The  inversion  for  proceeds  in  the  same  manner  for  |0|  <6C,  but  for  \8\  >  8C  the 
Cagniard  contour  must  be  deformed  around  a  branch  cut  on  the  real  axis  (cf.  [14], 
Fig.  7.12).  The  angular  variation  of  the  discontinuity  in  t/>,u  across  the  circular  part 
of  the  S-wavefront  is  then  obtained  by  evaluating  Re  [s2j3*F]  at  the  point  where  the 
Cagniard  contour  meets  the  branch  cut  in  the  upper  half-plane  (0  denotes  a  branch  of 
(62—  x2)!  with  positive  real  part).  Also,  the  part  of  the  inversion  contour  around  the 
branch  cut  gives  rise  to  a  head-wave,  which  precedes  the  5-wavefront,  and  it  con  tri¬ 
tributes  a  logarithmic  discontinuity  given  by 


(I—  v)Ko  v 
Ann  (nr)1  I  —  bv  cos  8 


L(8)  log  \t  —  br  1 , 


L(8)  =  {sin  28  ( a/b  +  cos  0)*/(c  +  b  cos  0)}  1m  [1/S  (b  cos  0)] . 


(A. 10) 


(A. 11) 


(vi)  The  stress  discontinuities  can  now  be  derived  from  the  relations  between  the  stress 
components  and  the  displacement  potentials  ([14],  §3.13),  by  using  the  operator 
identity  5/i )r  =  — (1  /c)  D/< >/,  which  holds  at  the  wavefronts  ([14],  §4.5.3;  c  =  cl,  ct  at 
the  P  and  5-wavefronts  respectively).  Thus,  for  the  normal  stress  a„  at  the  P- wave  front, 
we  have 

arr  =  (A  +  2/u.)  a2<f>,tt  -  (A a/r)  <f>,t  +  r  2</>,et , 

since  <fi  —0  before  the  5-wavefront  arrives.  In  view  of  the  wavefront  expansion  given 
for  <f>,tt  in  (A.9),  it  is  clear  that  the  wavefront  behaviour  of  o„  is  dominated  by  the 
step-discontinuity  in  the  term  (A  -f  2 p.)  a 2  </>,u.  This  leads  to  the  result  given  in  equations 
(1)  and  (2)  of  the  main  text.  The  behaviour  of  the  other  stress  components  is  deduced 
by  similar  considerations. 

3.  Displacement  of  the  crack  faces 
From  (A.l)  and  (A. 5)  we  have 

U„(x,  o,  s )  =  s-v*  CM-(x)l(d  +  x),  (A.  12) 

with  C  =  kKo/\/2 ■  The  inversion  with  respect  to  sy  is  easily  carried  out  for  x  >  0:  the  inversion 
contour  can  be  closed  to  the  left,  and  the  value  of  the  integral  is  simply  the  residue  at  x  =  —d. 
Thus, 

uy(x,  o,  s)  =  s  3'2  CM~(—d)  e~,dx, 

and  hence,  by  using  the  convolution  theorem  for  the  inversion  over  s,  we  find 

uy(x  >  o,  o,  t)  =  2ir_1  CM+(d)(t  —  xdy,  xd  <  t,  (A.  13) 

sine  M~(—d)  —  M+(d). 

The  inversion  for  x  <0  is  more  complicated:  the  inversion  contour  is  deformed  to  the  right 
onto  the  branch  cut,  and  the  resulting  integral  has  to  be  evaluated  numerically.  This  is  not  easy 
because  of  the  term  5  in  the  kernel  (cf.  (A. 3),  (A.4)).  One  can  use  a  rational  approximation  for 
5  (e.g.  [13],  Appendix),  but  we  have  preferred  instead  to  use  the  approximation  introduced 
in  1 ;  this  leads  to  an  inversion  in  closed  form.  The  idea  is  to  replace  (A. 3)  by 

A/+(x)  =  (h  +  x)!/(c  +  x)-  (A.  14) 

The  best  choice  for  h  is  discussed  in  1:  for  practical  purposes  one  can  take  h  =  0-956.  When 
the  inversion  contour  is  deformed  on  to  the  branch  cut  we  now  obtain 

u/x  <  o,  o ,  s)  =  s~3  2  (C>)  —  - -  e*l*  d£  . 

h  (d  +  ?)(?-  c) 

The  inversion  over  x  leads  to 

Uy(x  <  o,  o,  t)  =  2C(|jri/»r)!  F , 

where  F  is  a  convolution  integral  which  can  be  evaluated  in  closed  form  by  turning  it  into  a 
contour  integral  around  a  closed  loop  L  surrounding  a  branch  cut  between  h  !jc|  and  t  : 

1  J*  {(f  -  a  (6  Is!  -  w. 

2 nij,.a  +  \xd)a  -  cun)  4' 

If  the  loop  is  deformed  into  a  circle  about  the  origin  of  radius  r  -+oo,  the  value  of  the  integral 
can  be  derived  from  the  residues  at  cijr|,  ~\x\d,  oo,  when  t  <  r|x|.  For  t  >  c  |x|  the  pole  at 
C  =  clxt  is  not  crossed  in  the  process  of  deforming  the  contour  of  integration:  the  residues 
at  —1*1  d  and  oo  then  lead  to  the  Cauchy  principal  value  of  F.  We  thus  find,  with 

Rr  =  (c  -  6)!  (c\x\  -  r)V{(c  +  d)  |jri>) , 

Ri  =  (h  f  d) i  (t  +  \x\dY-H(c  +  d)  |jr|*| , 

Rr>  —  —  I  , 


that 


«»(—'/<■  <  x  <  o,o,  i)  =  2 C  (|jt| /n)1  ( Rc  +  Rd  +  R<*>) ,  (A.  15) 

Uy(—t/h  <  x  <  —lie,  o,  t)  =  2 C  (M/w)1  (/?<«  +  Ra>) .  (A. 16) 

Now,  to  determine  the  actual  displacement  of  the  crack  faces,  we  must  add  the  displacement 
given  by  (A.  13,  A.  15,  A.  16)  to  the  prior  static  crack  opening,  given  by  (see  [7]) 

UV(-I  <  X  <  O,  O,  I  <  0)  =  {o(l  -  v)/p}  {-x  (/  +  x)}1 . 

(Note  that  2 C  —  o(nl  )!  (1  —  v)//*)- 


4.  The  diffraction  of  the  starling  phase 

To  study  this  problem  it  is  more  convenient  to  suppose  that  it  is  the  crack  tip  at  x  =  0 
which  remains  stationary,  and  the  one  at  x  =  —I  which  moves  to  the  left  with  speed  v.  From  I : 
(30)  we  have 


f°° 

AK(s)-—  V2  oylf-u,s)di*(u,s)du. 

J  0 

Using  the  formula  for  the  Laplace  transform  of  a  product,  and  equation  (A. 5)  for  £yy  (allowing 
for  a  shift  of  origin),  we  obtain 


5  (AK/Ko)  = 


2  ni 


J  B(d  +  x)M< 


(X) 


dx. 


The  inversion  contour  B  is  deformed  to  the  right  on  to  the  branch  cut  to  give 


s  (AK/Ko)  = 


(d+O  MHO 


e  "t  d£, 


and  the  inversion  over  s  then  leads  to 


AK(n t  al)  = 


Ko 


Im  [M -(«//)] 


n\  a/(u  +  Id)  MHull) 


du  . 


(A.  17) 


The  essential  feature  is  a  pole  in  the  integrand;  it  is  more  clearly  evident  if  we  use  the 
approximate  kernel  (A.  14),  which  gives 


Ko  f'  /«+  le\  (u  —  hl\  du 


For  the  corresponding  anti-plane  problem  we  can  similarly  derive  the  result 

Ko 


AK(t>  bl)  = 


(A.  18) 


(A.  19) 
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Fig.  2 :  The  displacement  of  the  lower  half  of  the  crack  surface, 
y  ;  "  u»  y  1  ~  l/cL).  The  dashed  curve  shows  the 
initial  crack  opening.  A'  is  the  current  position  of  tip  A.  The 
depth  of  the  crack  relative  to  its  length  was  determined  by 
setting  o(1 -v)/2n  =  1. 


Fig.  3 :  Co-ordinates  and  notation  for  the  radiation  from  a  circular 
crack,  r  =  (x*-fy*)1'*. 


Fig.  4 :  The  step  response  A(t)  for  the  model  of  a  piezo-electric 
transducer,  with  W0/W,  =  2,  W0/W2  =  9. 


Fig.  5:  (a)  The  specimen  considered  in  Sectioi.  6,  showing  the 
location  of  transducers  A  and  6.  (b)  A  schematic  representation 
of  the  normal  surface  displacement  along  a  radial  line 
(8  =  30°)  on  z  =  h,  when  a  narrow  rectangular  P-pulse  is 
radiated  by  the  crack. 
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